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Our current understanding of the structure of nuclei with A < 8, including energy 
spectra, electromagnetic form factors, and weak transitions, is reviewed within 
the context of a realistic approach to nuclear dynamics based on two- and three- 
nucleon interactions and associated electro-weak currents. Low-energy radiative 
and weak capture reactions of astrophysical relevance involving these light systems 
are also discussed. 



1 Introduction 

Few-nucleon systems provide a unique opportunity for testing the simple, 
traditional picture of the nucleus as a system of point-like nucleons interact- 
ing among themselves via effective many-body potentials, and with external 
electro-weak probes via effective many-body currents. Through advances in 
computational techniques and facilities, the last few years have witnessed 
dramatic progress in numerically exact studies of the structure and dynamics 
of systems with mass number A < 8, including energy spectra of low-lying 
states, momentum distributions and cluster amplitudes, elastic and inelastic 
electromagnetic form factors, /3-decays, radiative and weak capture reactions 
at low energies, inclusive response to hadronic and electro-weak probes at 
intermediate energies. 

In the present talk, I will review the "nuclear standard model" outlined 
above, and present the extent to which it is successful in predicting some of the 
nuclear properties alluded to earlier. Of course, given the limited time, some 
of the theoretical and experimental developments will be treated cursorily. 
Nevertheless, I still hope to be able to convey a broad view of the intriguing 
and important studies in few-nucleon physics today. 



2 Potentials and Energy Spectra 

The Hamiltonian in the nuclear standard model is written as 
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i i<j i<j<k 

where the kinetic energy operator Kt has charge-independent and charge- 
symmetry-breaking components due to the difference in proton and neutron 
masses, and Vij and Vijk are two- and three-nucleon potentials. 

The two-nucleon potential consists of a long range part due to pion ex- 
change, and a short-range part parameterized either in terms of heavy meson 
exchanges as, for example, in the Bonn potential El, or via suitable operators 
and strength functions, as in the Argonne Vig (AV18) potential 13. The short- 
range terms in these potentials are then constrained to fit pp and np scattering 
data up to energies of ~ 350 MeV in the laboratory, and the deuteron bind- 
ing energy.. The modern models mentioned above provide fits to the Nijmegen 
data-base □ characterized by \ 2 P er datum very close to one, and should there- 
fore be viewed as phase-equivalent. The AV18 model is most widely used; it 
has the form 



= o& . ( 2 ) 

p=l,18 

where the first fourteen operators are isoscalar, 

0£ =1 - 14 = [1,£T< • <Tj,Sij, (L • S^I/SLV, • (TJ, (L • S)?-] ® [l,Ti • Tj] , (3) 

while the last four isospin-symmetry-breaking operators have isovector and 
isotensor character, 

Off 15 18 = Tij , crj ■ (TjTij , SijTij , (n + Tj) z . (4) 

Here Sij is the tensor operator, and is defined as = 3ri Z Tj Z — n ■ Tj. 
Unique among the modern potentials, the AV18 includes a fairly complete 
treatment of the electromagnetic interaction, since it retains, in addition to 
the leading Coulomb term, also contributions from magnetic moment inter- 
actions, vacuum polarization and two-photon exchange corrections. These 
terms, while typically very small (for example, in the deuteron the magnetic 
dipole-dipole interaction gives 18 keV extra repulsion □), need to be taken 
into account when very accurate predictions are required, as in the case, for 
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example, of studies of energy differences of isomultiplet states □, or the cross 
section for proton weak capture on proton at keV energies u. 

It is now well established that two-nucleon potentials alone underbind 
nuclei □: for example, the AV18 and Bonn models give Q, in numerically exact 
calculations, binding energies of 24.28 MeV and 26.26 MeV respectively, which 
should be compared to the experimental value of 28.3 MeV. Moreover, 6 Li and 
7 Li are unstable against breakup into ad and at clusters, respectively, and that 
energy differences are not, in general, well predicted, when only two-nucleon 
potentials are retained in the Hamiltonian. 

Important components of the three-nucleon potential arise from the in- 
ternal structure of the nucleon. Since all degrees of freedom other than the 
nucleon have been integrated out, the presence of virtual A resonances, for 
example, induces three-nucleon potentials. They are written as 



ViA 



vijk = V$ + V$ H , (5) 
where V 27T is the "long-range" term, resulting from the intermediate excitation 
of a A with pion exchanges involving the other two nucleons, known as the 
Fujita-Miyazawa term 13. This term is present in all models, such as the 
Tucson-Melbourne potential □ or the series of Urbana models The Urbana 
models parameterize V R as 



V R = A R J- T%(rij)T%(rjk) , (6) 

cyclic ijk 

where T^ir) is the strength function of the pion-exchange tensor interaction. 
This term is meant to simulate the dispersive effects that are required when 
integrating out A degrees of freedom. The strengths of the Fujita-Miyazawa 
and dispersive terms are then determined, in the Urbana models, by fitting 
the triton binding energy and the saturation density of nuclear matter. 

The Hamiltonian consisting of the AV18 two-nucleon and Urbana-IX 
three-nucleon potentials (AV18/UIX) predicts reasonably well the low-lying 
energy spectra of systems.. with A < 8 nucleons in "exact" Green's function 
Monte Carlo calculations □. The experimental binding energies of the a par- 
ticle is exactly reproduced, while those of the A=6-8 systems are underpre- 
dicted by a few percent. This undcrbinding becomes (relatively) more and 
more severe as the neutron-proton asymmetry increases. An additional failure 
of this Hamiltonian model is the underprediction of spin-orbit splittings in the 
excitation spectra of these light systems. These failures-have in fact led to the 
development of new three-nucleon interaction models E2I. These newly devel- 
oped models, denoted as Illinois models, incorporate the Fujita-Miyazawa and 
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dispersive terms discussed above, but include in addition multipion exchange 
terms involving excitation of one or two A's, so-called pion-ring diagrams, as 
well as the terms arising from S-wave pion rescattering, required by chiral 
symmetry. 

3 The Nuclear Electromagnetic Current 

The nuclear current operator consists of one- and many-body terms that op- 
erate on the nucleon degrees of freedom: 

j(q) = £j« (1) (q)+£jg ) (q) + £ J$(q). (7) 

i i<j i<j<k 

where q is the momentum transfer, and the one-body operator has the 
standard expression in terms of single-nucleon convection and magnetiza- 
tion currents. The two-body current operator has "modeLindependent" and 
"model-dependent" components (for a review, see Ref. B). The model- 
independent terms are obtained from the charge-independent part of the 
AV18, and by construction satisfy current conservation with this interac- 
tion. The leading operator is the isovector "7r-like" current obtained from the 
isospin-dependent spin-spin and tensor interactions. The latter also generate 
an isovector "p-like "current, while additional model- independent isoscalar 
and isovector currents arise from the central and momentum-dependent inter- 
actions. These currents are short-ranged and numerically far less important 
than the 7r-like .current. Finally, models for three-body currents have been 
derived in Ref. £3, however the associated contributions have been found to 
be very small in studies of the magnetic structure of the trinucleons L3. 

The model-dependent currents are purely transverse and therefore can- 
not be directly linked to the underlying two-nucleon interaction. Among 
them, those associated with the A-isobar are the most important ones in the 
momentum-transfer regime being discussed here. These cuxrents are treated 
within the transition-correlation-operator (TCO) scheme a scaled-down 
approach to a full A+A coupled-channel treatment. In the TCO scheme, the 
A degrees of freedom are explicitly included in the nuclear wave functions by 
writing 



N+A 



(8) 



bologna 'new: submitted to World Scientific on February 8, 2008 4 



where ^ is the purely nucleonic component, S is the symmetrizer and the 
transition correlations U^ R are short-range operators, that convert NN pairs 
into iVA and A A pairs. In the results reported here, the \& is taken from CHH 
solutions of the AV18/UIX Hamiltonian with nucleons only interactions, while 
the Ufj R is obtained from two-body bound and low-energy scattering state 
solutions of the full N-A coupled-channel problem. Both jNA and 7AA Mi 
couplings are considered with their values, n-yNA = 3 n.m. and /i 7 AA = 4.35 
n.m., obtained from data Ej. 



4 The pd Radiative Capture 



There are now available many high-quality data, including differential cross 
sections, vector and tensor analyzing powers, and photon polarization coef- 
firipnts nri tlip pd radiative capture at cm. energies ranging from to 2 
MeV EaliSESIiZl. These data indicate that the reaction proceeds predomi- 
nantly through S- and P-wave capture. The aim here is to verify the extent 
to which they can be described satisfactorily by a calculation based on a real- 
istic Hamiltonian (the AV18/UIX model) and a current operator constructed 
consistently with the two- and three- nucleon interactions £3. 

The predicted angular distributions of the differential cross section a (9), 
vector and tensor analyzing powers A y {9) and T 2 q{9), and photon linear po- 
larization coefficient P-y(9) are compared with the TUNL data below 50 keV 
from Refs. E_3'E_3 in Fig [l| The agreement between the full theory, including 
many-body current contributions, and experiment is generally good. However, 
a closer inspection of the figure reveals the presence of significative discrepan- 
cies between theory and experiment in the smalLangle behavior of u{9) and 
22o(0), as well as in the 5-factor below 40 keV The S-wave capture pro- 
ceeds mostly through the M\ transitions connecting the doublet and quartet 
pd states to 3 He-the associated reduced matrix elements (RMEs) are denoted 
by m,2 and 7714, respectively. The situation for P-wave capture is more com- 
plex, although at energies below 50 keV it is dominated by the E\ transitions 
from the doublet and quartet pd states having channel spin 5=1/2, whose 
RMEs I denote as p 2 and p^. The E\ transitions involving the channel spin 
5 = 3/2 states, Bihile smaller, do play an important role in T-^iG). 

The TUNL M and Wisconsin 113 groups have determined the leading Mi 
and Ei RMEs via fits to the measured observables. The results of this fitting 
procedure are compared with the calculated RMEs in Table |l|. The phase of 
each RME is simply related to the elastic pd phase shift LLJ, which at these low 
energies is essentially the Coulomb phase shift. As can be seen from Table |l|, 
the most significant differences between theoretical and experimental RMEs 
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d+p Capture 
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Figure 1. The energy integrated cross section cr(9)/ao (47rao is the total cross section), 
vector analyzing power A y (S), tensor analyzing power T2o(d) and photon linear polariza- 
tion coefficient Pj(9) obtained with the AV18/UIX Hamiltonian model and one-body only 
(dashed line) or both one- rmd many-body (solid line) currents are compared with the 
experimental results of Rcf. EJ. 



are found for \p4,\. The theoretical overprediction of p4 is the cause of the 
discrepancies mentioned above in the low-energy (< 50 keV) S- factor and 
small angle a (6). 

It is interesting to analyze the ratio tei = \pa/p2\ 2 - Theory gives tei — 1, 
while from the fit it results that tex ~ 0.74 ± 0.04. It is important to stress 
that the calculation of these RMEs is not influenced by uncertainties in the 
two-body currents, since their values are entirely given by the long-wavelength 
form of the E\ operator (Siegert's theorem), which has no spin-dependence 
(for a thorough discussion of the validity of the long- wavelength approximation 
in Ei transitions, particularly suppressed ones, see Ref. ta). It is therefore of 
interest to examine more closely the origin of the above discrepancy. If the 
interactions between the p and d clusters are switched off, the relation vex — 1 
then simply follows from angular momentum algebra. Deviations of this ratio 
from one are therefore to be ascribed to differences induced by the interactions 
in the 5=1/2 doublet and quartet wave functions. The AV18/UIX interactions 
in these channels do not change the ratio above significantly. It should be 
emphasized that the studies carried out up until now ignore, in the continuum 
states, the effects arising from electromagnetic interactions beyond the static 
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Coulomb interaction between protons. It is not clear whether the inclusion of 
these long-range interactions, in particular their spin-orbit component, could 
explain the splitting between the P2 and 774 RMEs-pbserved at very low energy. 
This discrepancy seems to disappear at 2 MeV 113. 

Table 1. Magnitudes of the leading Mi and E\ RMEs for pd capture at E p = 40 keV. 



RME IA FULL FIT 



\m 2 \ 


0.172 


0.322 


0.340±0.010 


|m 4 | 


0.174 


0.157 


0.157±0.007 




0.346 


0.371 


0.363±0.014 


M 


0.343 


0.378 


0.312±0.009 



Finally, the doublet 777,2 RME is underpredicted by theory at the 5 % 
level. On the other hand, the cross section for nd capture at thermal neutron 
energy is calculated to be 578 fib with the AVLS7UIX model, which is 15 % 
larger than the experimental value (508±15) fib ta. Of course, Mi transitions, 
particularly doublet ones, are significantly influenced by many-body current 
contributions. Indeed, an analysis of_the isoscalar (/is) an d isovector (juy) 
magnetic moments of the trinucleons 113 suggests that the present model for 
the isoscalar two-body currents, constructed from the AV18 spin-orbit and 
quadratic-momentum dependent interactions, tends to overestimate /is by 
about 5 %. The experimental value for /jy, however, is almost perfectly 
reproduced. It appears that the present model for two-body currents needs 
to be improved. 

5 The Nuclear Weak Current and the j7 3 He Weak Capture 

The nuclear weak current and charge operators consist of vector and axial- 
vector parts, with corresponding one- and many-body components. The 
weak vector current and charge are constructed from the corresponding 
(isovector) electromagnetic terms, in accordance with the conserved-vector- 
current hypothesis, and thus have cJ "model-independent" and "model- 
dependent" components. The former are determined by the interactions, the 
latter include the transverse currents associated with A excitation. 

The leading many-body terms in the axial current, in contrast to the 
case of the weak vector (or electromagnetic) current, are those due to A 
excitation, which are treated within the TCO scheme, discussed above. The 
axial charge operator includes the long-range pion-exchange term Ell, required 
by low-energy theorems and the partially-conserved-axial-current relation, as 
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well as the (expected) leading short-range terms constructed fmin the central 
and spin-orbit components of the nucleon-nucleon interaction c3. 

The largest model dependence is in the weak axial current. The iVA axial 
coupling constant g* A is not well known. In the quark-model, it is related to 
the axial coupling constant of the nucleon by the relations g* A = (6\^2/5)gA- 
This value has often been used in the literature in the calculation of A-induced 
axial current contributions to weak transitions. However, given the uncertain- 
ties inherent to quark-model predictions, a more reliable estimate for g* A is 
obtainecLby determining its value phenomenologically. It is well established 
by now u that one-body axial currents lead to a ~ 4 % underprediction of 
the measured Gamow- Teller matrix element in tritium /3-decay. This small 
4 % discrepancy can then be used to determine g* A Ej. While this procedure 
is inherently model dependent, its actual model dependence is in fact very 
weak, as has been shown in Ref. 

The calculated values for the^strophysical S-factor in the energy range 
0-10 keV are listed in Table || E3. Inspection of the table shows that: (i) 
the energy dependence is rather weak, the value at 10 keV is only about 4 
% larger than that at keV; (ii) the P-wave capture states are found to be 
important, contributing about 40 % of the calculated S'-factor. However, the 
contributions from D-wave channels are expected to be very small, as explicitly 
verified in 3 D! capture, (iii) The many-body axial currents associated with A 
excitation play a crucial role in the (dominant) 3 Si capture, where they reduce 
the S'-factor by more than a factor of four; thus the destructive interference 
between the one- and many-body current contributions, obtained in Ref. is 
confirmed in the study of Ref. E3, based on more accurate wave functions. The 
(suppressed) one-body contribution comes mostly from transitions involving 
the D-state components of the 3 He and 4 He wave functions, while the many- 
body contributions are predominantly due to transitions connecting the S- 
state in 3 He to the D-state in 4 He, or viceversa. 



Table 2. The hep S'-factor, in units of 10 -20 keV b, calculated with CHH wave functions 
corresponding to the AV18/UIX Hamiltonian model, at p 3 He cm. energies E=0, 5, and 
10 keV. The rows labelled "one-body" and "full" list the contributions obtained by retaining 
the one-body only and both one- and many-body terms in the nuclear weak current. The 
contributions due the 3 Si channel only and all S- and P-wave channels are listed separately. 





E=0 keV 


E=5 keV 


E=W keV 




3 Si S+P 


3 Si S+P 




S+P 


one-body 
full 


26.4 29.0 
6.38 9.64 


25.9 28.7 
6.20 9.70 


26.2 
6.36 


29.3 
10.1 
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ncement, while very significant, is smaller than that first suggested 
J j ' Even though this result is inherently model dependent, it is un- 



The chief conclusion of Ref. El is that the hep 5-factor is predicted to be^p 
4.5 times larger than the value adopted in the standard solar model (SSM) e3. 
This en 
in Ref. 

likely that the model dependence is large enough to accommodate a drastic 
increase in the value obtained here. Indeed, calculations using Hamiltonians 
based on the AV18 two-nucleon interaction only and the older AV14/UVIII 
two- and three-nucleon interactions cB predict zero energy S-factor values of 
12.1 x 1CT 20 keV b and 10.2 x 1CT 20 keV b, respectively. It should be stressed, 
however, that the AV18 model, in contrast to the AV14/UVIII, does not repro- 
duce the experimental binding energies and low-energy scattering parameters 
of the three- and four-nucleon systems. The AV14/UVIII prediction is only 
6 % larger than the AV18/UIX zero-energy result. This 6 % variation should 
provide a fairly realistic estimate of the theoretical uncertainty due to the 
model dependence. The precise calculation of the ^-factor and the conse- 
quent absolute prediction for the hep neutrino flux should allow much greater 
discrimination among proposed solar neutrino oscillation solutions Ej. 



6 Conclusions and Acknowledgments 

Improvements in the modeling of two- and three-nucleon interactions and 
nuclear electro- weak currents, and the significant progress made in the last few 
years in the description of bound and continuum wave functions, make it now 
possible to perform first-principle calculations of interesting nuclear properties 
of light nuclei. Experimentally known electromagnetic and weak transitions of 
systems in the mass range 2 < A < 8 provide powerful constraints on models 
of nuclear currents. 
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pande, S.C. Pieper, D.O. Riska, S. Rosati, M. Viviani, and R.B. Wiringa for 
their many important contributions to the work reported here. This work was 
supported by DOE contract DE-AC05-84ER40150 under which the Southeast- 
ern Universities Research Association (SURA) operates the Thomas Jefferson 
National Accelerator Facility. 
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